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ABSTRACT 



The Hyades cluster is known to harbour ten so-called classical white dwarf members. Numerous studies through the years have 
predicted that more than twice this amount of degenerate stars should be associated with the cluster. Using the PPMXL catalog of 
proper motions and positions, a recent study proposed 17 new white dwarf candidates. We review the membership of these candidates 
by using published spectroscopic and photometric observations, as well as by simulating the contamination from field white dwarfs. 
In addition to the ten classical Hyades white dwarfs, we find six white dwarfs that may be of Hyades origin and three more objects 
that have an uncertain membership status due to their unknown or imprecise atmospheric parameters. Among those, two to three are 
expected as field stars contamination. Accurate radial velocity measurements will confirm or reject the candidates. One consequence 
is that the longstanding problem that no white dwarf older than ~340 Myr appears to be associated with the cluster remains unsolved. 

Key words, open clusters and associations: individual: Hyades — stars: atmospheres — white dwarfs 



1. Introduction 

The Hy ades cluster, with a center-of-mass distance of 46.45 + 
0.5 pc (Ivan Leeuwenll2009t) and an estimated age of 625 + 50 
Myr, derived from the heliu m abundance and isochrone mod- 
eling (Perrvm an et al.lll998l), is certainly one o f the most stud- 
ied open clusters. Recently, Ros er et all d201 ll hereafter Paper 
I) identified 724 probable main-sequence members within a ra- 
dius of 30 pc from the cluster center. Their samp le is based 
on a k inematic selection from the PPMXL Catalog (Rose r et"a"D 
2010) with the convergent point method, further constrained by 
photometric data. They determined individual masses for their 
candidates and, among others, studied the mass distribution as a 
function of the distance to the center. It was confirmed that the 
present day tidal radius of the cluster is ~9 pc, and that this value 
would only change by a few percent under different assumptions 
on the fraction of binaries, a quantity that is currently not well 
constrained. 

In this work, we restrict our study to Hyades white dwarfs, 
even for which the literature is extensive. It was already con- 
cluded long ago that old white dwarfs with an age greater 
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than ~300 Myr were missing from the cluster (lTinslevlll974t 
Ivan den Heuvell 1 19751). and that th e y might have eva p orated 
dPels et al.l 1 1975HWeidemannl 119771) . IWeidemann et all d 19921) 
and lEggenl (1 19931) further constrained this scenario by suggest- 
ing that old evaporated Hyades white dwarfs form long tidal 
tails for which the density is much lower than that of field 
stars. To this day, it is still the case that only 10 constituents 
are confirmed, and N-body simulations systematically predict 
that at l east two times more white dwarfs were form ed in the 
Hyades (Porte gies Zwart et alfeOOlUErnst et~ai1l201 ll) . The idea 



that many remnants evaporated from the cluster remains a likely 
scenario, which is further suggested by the fact that Paper I 
found prominent tidal tails for main-sequence stars, with almost 
as many members currently leaving the cluster than members 
bound to the cluster. However, it is also possible that the miss- 
ing white dwarfs were simp ly never formed beca use of a very 
steep initial mass function ( Boh m- Vi tense 20071) o r that they 
are hi dden in binaries (Portegies Zwa rt et al. 2001). IWilliamsl 
(2004) tried to bring the observed number of white dwarfs in 
the Hyades into agreement with calculated numbers by vary- 
ing the slope of the IMF, the binary fraction, the binary mass 
ratio and the effect of dynamical evolution. He concludes that, 
from simulations alone, it is not possible to differentiate between 
the above scenarios. A deficiency of white dwarfs has also been 
observed in several other o pen clusters (see, e.g.. iRicher et al.l 
119981: iKahrai et al.l l2001bt Williams & Boltd 120071) . although 
there are also many instances where the observed and predicted 
numbers of degenerate objects are in agreement ( Kalir ai et al.l 
1200 lal 120031: iRubin et al.ll2008l) . 

In a recent analysis, Schilb ach & RoseJ d2012l hereafter 
Paper II) extended the systematic search of Hyades members to 
white dwarfs up to 40 pc from the cluster center. The authors did 
recover all 10 classical components, and further identified 17 
white dwarfs as new candidates, all of them outside of the tidal 
radius. If these degenerate stars are confirmed as members, they 
could potentially explain the discrepancy between the observed 
and predicted number of associated objects, and also confirm 
that white dwarfs are slowly leaving the cluster. 

The 10 classical Hyades white dwarfs are certainly some 
of the most well known and observed hot remnants. Given 
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that new atmospheric parameters have recently been published 
for most of these ob j ects based on spectr oscopic analyses 
dGianninas et al.l 1201 U iBergeron et al.l 1201 If) reiving o n im- 
proved model atmospheres dTremblav & Bergeron 2009), it is 
certainly pertinent to review their properties as Hyades represen- 
tatives. Furthermore, most of the 17 new candidates identified 
in Paper II also have well determined atmospheric parameters, 
from spectroscopic or photometric analyses. Consequently, the 
next step to constrain their membership is to compare quantities 
derived from kinematic studies, e.g. kinematic distances, with 
those derived from their atmospheric properties, which includes 
the so-cal led spectroscopic d istance, based on predicted lumi- 
nosities dHolberg et al.l 12008). The spectroscopic distances to- 
gether with the proper motions can be used to update the space 
velocities for candidates, i.e. the residual 2D velocity with re- 
spect to the Hyades bulk motion, which is important to under- 
stand the structure of the cluster. 

This work is structured as follow. In Sect. 2, we describe the 
Hyades white dwarfs sample of Paper II and Sect. 3 follows with 
a description of analyses of their spectra and photometry. Based 
on their derived atmospheric parameters, we compute, in Sect. 
4, spectroscopic distances for all aspirant Hyades remnants. In 
Sect. 5, we derive improved residual velocities, and try to un- 
derstand the global status of new candidates by simulating the 
field white dwarfs contamination. Then, Sect. 6 continues with a 
discussion on the membership of new candidates and the impli- 
cations on the Hyades evolution. The summary is presented in 
Sect. 7. 



2. Kinematic and photometric Hyades candidates 

The starting point of our analysis are the Hyades white dwarf 
candidates identified in Paper II. The selection process is ex- 
plained in details in Paper I and II. In brief, the primary 
data of observations are the proper motions and positions in 
the PPMXL catalog. The PPMXL is an all-sky survey with 
a limiting magnitude of V ~ 20 including about 900 mil- 
lion objects. The PPM XL data was cr oss-matched and com- 
bined with the UCAC3 dZacharias et alJl2010l) and the CMC 14 
dCopenhagen Univ. Obs. et al. I I2006I) . and referred as a whole as 
the Carlsberg-UCAC (CU) subset. Along with the kinematics, 
this data set also includes 2MASS JHK S and CMC 14 r' magni- 
tudes. 

Candidates are first selected using the convergent point 
method (Paper I). More precisely, in the case of white dwarf 
candidates, an upper limit of 5 km s _1 is allowed for the dif- 
ference between their tangential motion and the bulk motion 
of the Hyades. The convergent point method then predicts, for 
each possible member, a radial velocity and a distance, hereafter 
named the kinematic distance (akin), based only on the position 
(a, 6) and proper motion (ju a , fig). To avoid contamination by 
field objects, only the new candidates with a ifkj n prediction that 
put them closer than 40 pc from the cluster center are kept. 

The second step used in Paper II to construct the Hyades 
sample was to cross-match all kinematically suited candidates 
with white dwarf catalogues, suc h as the online version of the 
Villanova White Dwarf Catalosfl (McC ook & Sionll 19991) . From 
this procedure, the authors identified a total of 37 white dwarfs 
as kinematic candidates. We note that at this point, a fair num- 
ber of these objects could be background objects with proper 
motions that mimic the Hyades motion. Therefore, the final step 
was to examine the available photometric data of the candidates. 



These white dwarfs are unfortunately at the detection limit of 
the 2MASS and CMC 14 survey. Even if JHK S and r' mag- 
nitudes are available for most candidates, a high scatter was 
found in the J - Ks colors, which is consistent with the conclu- 
sion of lTrembla v & Bergeron! d2007l) that one should be cautious 
with K s magnitudes near the faint detection limit of the survey. 
Therefore, it was opted to include in the analysis B and V mag- 
nitudes from the literature even if the accuracy of the available 
BV was far from being uniform and varied from several mil- 
limag (photoelectric measurements) to several tenths of a mag 
(photographic measurements). 

Using the kinematic distance and an observed magnitude, 
a "kinematic" absolute magnitude (kin. M) can be computed. 
In Fig. Q] we reproduce the kin. My vs. B - V and kin. Mj vs. 
r'—J diagrams made in Paper II, but now compared with theoret- 
ical w hite dwarf cooling sequences (Wood 1995; Fontaine et al.l 
2001). True Hyades members are expected to follow the theo- 
retical sequences. On the contrary, background and foreground 
stars are expected to be shifted vertically in relation to the the- 
oretical sequences, because their estimated o^in differs from the 
true distance. We can see from Fig. Q] that the 10 white dwarfs 
identified as probable non-members in Paper II (open circles) 
lie below the theoretical cooling sequences. They would have to 
be super massive (log g > 9.0) to be cluster members. Instead, 
they are much more likely distant objects mimicking the Hyades 
proper motion. 
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Fig. 1: Kinematic absolute magnitude versus color diagrams for 
Hyades white dwarf candidates (filled circles) and probable non- 
Hyades members (open circles) identified in Paper II. Binaries 
are marked as filled triangles and by their numbers from Table 1, 
although in most cases, they are too bright for the range of ab- 
solute magnitudes shown in the figure. Theoretical cooling se- 
quences are shown from log g = 7.0 (dotted line) to log g = 9.0 
(solid line) in steps of 0.5 dex. 



All 10 of the so-called classical Hyades members were re- 
covered in Paper II. Removing 10 white dwarfs identified as 
probable non-members, as discussed above, there remain 17 
new candidates, all lying outside of the tidal radius of ~9 pc. 
However, five of the remaining candidates are more than 20 
pc away from the cluster center in the Z direction (perpendic- 
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ular to the galactic plane) and main-sequence stars with known 
radial velocities in these regions have been ruled out as mem- 
bers (Paper I). We note that three of the classical members and 
two of the new candidates are unresolved binaries with a main- 
sequence companion. 

2.1. New selection of Hyades candidates 

The kinematic distances of the Hyades white dwarf candidates 
were obtained in Paper II by the requirement of the minimum 
difference between their tangential velocity and the cluster mo- 
tion, with an upper limit of the residual velocity of 5 km s _1 . 
Later in this work (see Sect. 5 below) we recalculate residual 
velocities of the candidates by assuming the spectroscopic dis- 
tances derived in Sect. 4. We found that most of the candidates 
outside of the tidal radius have residual velocities up to 8 km s _1 . 
Therefore, we conclude that an upper limit of 5 km s for the 
residual velocity may have been a too rigorous requirement for 
detecting escaping white dwarfs. 

To obtain a complete sample of possible Hyades candidates, 
we revisited the selection process of Paper II by using a higher 
upper limit of 8 km s for the residual velocity. We found 14 ad- 
ditional white dwarfs for which the kinematic requirement was 
fulfilled. All but three candidates were later rejected by spectro- 
scopic criteria. The three promising candidates are listed as star 
Nos. 38 to 40 at the bottom of Table 1 . 



3. Spectroscopic and photometric stellar 
parameters 

An improved diagnostic of the Hyades membership of white 
dwarfs can be made based on known stellar parameters. In 
Table 1, we have compiled the r e ff, log g and total stellar masses 
found in the literature for 37 of the 40 candidates. The corre- 
sponding absolute V magnitudes (My) and white dwarf cool- 
ing time (log r) have then been computed using the evolution- 
ary models with thick layers (hydrogen-rich atmosphere DA 
white dwa rfs) and thin layers (he lium-rich atmosphere DB white 
dwarfs) of Fontaine et al. (2001, carbon-oxygen cores) for T e s < 
30,000 K and IWoodl d!995l carbon cores) above this tempera- 
ture. When available in the cited literature, we also compiled the 
distances, although we note that the calculation of improved dis- 
tances is the goal of our Sect. 4 below. Omitted from Table 1 
is the binary white dwarf WD 0217+375 (No. 13), which was 
designated as an M5 V+D A pair without a de tailed spectral anal- 
ysis of the white dwarf dSilvestri et al.l 120051) . Also omitted is 
LP 649-0071 (No. 12), a M dwarf star with a suspected but 
spectroscopically unconfirmed white dwarf companion , previ- 
ously identified in Luyten's White Dwarf Catalogues dLuvtenl 
Il977l) . Finally, the unconfirmed white dwarf and probable non- 
member 1RXSJ062052.2+132436 (No. 34) also has no precise 
atmospheric parameters. 

The majority of the candidates are hydrogen-rich DA white 
dwarfs. The most successful and precise technique to identify 
the atmospheric parameters of DA white dwarfs is to fit the 
observed Balmer line profiles with synthetic profiles predicted 
from m odel atmospheres, which is cal l ed the spectroscopic tech- 
nique dWeidem ann & K oesterl Il980t ISchulz & Wegnerl 1 198 it 
Berg eron et alj 1992b . As much as 26 objects in Table 1 rely on 
the spectroscopic analysis of the Balmer lines to identify their 
atmospheric parameters (references 2, 7, 11 and 12 in the table). 
Among them, 22 are f rom the White Dwarf Catalog survey of 
iGianninas et al.l d201 ll) and two are from the updated SDSS anal- 



ysis of Trembl ay et al.l (1201 lal) . Both of these analyses adopted 
the m ost recent model atmospheres of iTremblav & Bergeron! 
d2009l) . relying on Balmer line profiles including non-ideal gas 
effects. The quoted uncertainties in Table 1 are from the fitting 
procedure. However, more realistic estimates of the uncertainties 
can b e made by fi t ting di fferent observations of the same objects, 
and iLiebert et all d2005l) found spreads of 1.2% and 0.038 dex 
in r e ff and log g, respectively. For the remaining of this work, 
we use these values as the minimum uncertainties on the atmo- 
spheric parameters. We note that possible systematic uncertain- 
ties from the models are not included. Yet, promising Hyades 
candidates (Nos. 1-27) with spectroscopic parameters in Table 1 
have radiative atmosphere s and are outside the range of the so- 
called high-mass problem dKoester et al.ll20 09a; Trembla y et al.l 
1201 lbl) . For these objects (T e s > 13,000 K), the atmospheric 
parameters are in good agreement with the few observational 
co nstraints, such as parallax measurements (see, e.g., Fig. 14 
of ITremblav & Bergeronll2009l) . We emphasize that the current 
work can also serve as a validation for recent DA models since 
the distanc e to the Hyades cen ter is fairly well constrained via 
Hipparcos dvan Leeuwenll2009l) . Finally, the helium and hydro- 
gen lines in the DBA star WD 0437+1 38 have been fitted spec- 
troscopically by iBergeron et alj d201 lb with improved helium- 
rich models. 

Among the remaining candidates, WD 0433+270 (No. 21) 
is a cool DA with only a faint Ha line for which a spec- 
troscopic fit is not possible. The BVRIJHK photo metric en- 
ergy d istribution of this object has been fitted by Ber geron et al.l 
(2001) to determine T e ff and also the \ogg using the con- 
straint fro m the trigonometric paralla x. We use here the re- 
sults from iGiammichele et all d2012l) who refitted the data 
using models with the Ly man a H-H2 neutral broadening 
(Kowalski & Saumonll2006h . The best fit parameters were used 
to predict a Ha line profile that is in close agreement with the 
observed one (see Fig. 1 If from Giammichele et al.), and there- 
fore we can consider that the atmospheric parameters for this ob- 
ject are similarly precise as those found with the spectroscopic 
method for the previously discussed stars. The cool and nearby 
magnetic DQP WD 0548 -001 (No. 39) has also been fitted by 
IGiammichele et all (2012) using a similar method, although with 
an additional free parameter due to the presence of carbon in the 
helium-rich atmosphere. All in all, 29 of the 37 objects in Table 1 
are typical DA, DB or DQ with well constrained stellar param- 
eters. We remind the reader, however, that two magnetic objects 
(Nos. 24 and 39) and two objects that have been fitted without 
the improved Stark broadening profiles (Nos. 20 and 35) likely 
have underestimated uncertainties. 

Other objects in Table 1 have more uncertain atmospheric 
parameters. WD 0816+376 (No. 27) is a DAP white dwarf with 
a high magnetic field that is part of the Kiso survey sample of 
iLimo ges & Bergeron (2010). Even if the obj ect is include d in 
their final sample, they used the estimation of lJordarJ d 1993b for 
the temperature and fixed the gravity at \ogg = 8.0. Therefore, 
we chose a characteristic uncertainty of 0.5 dex for the grav- 
ity of this object. The situation is similar for the poorly stud- 
ied cool DC remnant WD 0648+368 (No. 40) and the DQA 
SDSS J0236 14.44-080804.9 (probable non-member) for which 
we also fixed the gravity. 

Of the 4 objects in Table 1 that are unresolved binaries (Nos. 
1, 4, 7 and 36), WD 0429+176 and WD 0816+387 can be fit- 
ted relatively well using the spectroscopic method by removing 
iteratively the M dwarf contribution. In contrast, the classical 
remnants WD 0347+171 and WD 0418+137 have F- and K- 
type main-sequence companions that completely dominate the 



3 



Tremblay et al.: Spectroscopic and photometric studies of white dwarfs in the Hyades 



optical flux. As a consequence, these two white dwarfs were 
observed in the UV to constrain their atmospheric parameters, 
and the spectroscopic method applied to the Lyman lines is not 
as well documented and precise as for the Balmer line analysis 
dVennes et al.ll2005h . For one of them, the gravity could not be 
constrained from the Lyman a spectrum alone, and we chose a 
characteristic value of log g — 8.0 + 0.5. 

We have at hand, from the PPMXL and paper II, a photo- 
metric data set with Johnson BV, CMC 14 r' and 2MASS JHK S 
magnitudes for most objects in Table 1 . We therefore proceeded 
to fit the photometric data with grids of pure-H and pure-He 
white dwarf models, with and the solid angle R/D, where 
R is the stellar radius and D the distance from the Sun, as free 
parameters. This method is not sensitive to the gravity, and we 
used a fixed value of loggf = 8.0. Consequently, this method 
can only be used to give an independent valu^3 of T e s, called 
photometric temperature be low. We rely on the prescription of 
Holberg & Bergeron (2006) to convert the magnitudes into aver- 
age fluxes. 

From our photometric fits, we find that in all but one case 
the photometric temperature agrees with the spectroscopic value. 
The exception is WD 0231-054 for which the JHK flux ap- 
pears too high and suggests a lower for this object. We 
have no clear explanation for this behaviour. Most objects in 
Table 1 are hot and blue, and their colors, in the Rayleigh-Jeans 
regime, are not very sensitive to (or the absolute magnitude, 
as can be seen in Fig. [TJ. Therefore, the photometric tempera- 
tures are less precise and we rely, whenever possible, on spec- 
troscopically derived temperatures. There are four cases where 
no spectroscopic parameters are available. The first is the DA 
WD 0120-0 24 (No. 11), which h as only been fitted using the 
V -I color in lSilvestri et al.l (1200 lb . On the contrary, our data set 
includes 2MASS JHKs photometry, and w e could also find for 
this object the Johnson BVRI photometry (Smith1 [T997h . Our fit 
is shown in F ig. [2] with the gravity fixed at the value found by 
ISilvestri et al.1 (120011) from a gravitational redshift measurement. 
We note a discrepancy at the K$ band, which is not unusual as 
it is close to the detection limit of 2MASS. The atmospheric pa- 
rameters should be considered as relatively well constrained for 
this object, despite the fact that masses determined from gravita- 
tional redshifts are generally not as preci se as those found fro m 
spectroscopic or parallax measurements dBergeron et al.l l2007). 

In Fig. [3] we show our BV and 2MASS JH fit of LP 246- 
14 (No. 15), which is found in the Luyten's White Dwarf 
Catalogues but is not a spectroscopically confirmed white dwarf, 
and no atmospheric parameters have yet been published. The fit 
is consistent with a relatively cool white dwarf, but the lack of a 
spectrum for this object prevents us to conclude that it is a white 
dwarf. While in the next section it is shown that the d^m value 
is consistent with a gravity of log g ~ 8, it does not exclude 
that it is a background lower gravity object with a high tangen- 
tial velocity. Therefore, the atmospheric parameters presented in 
Table 1, an average of the pure-H and pure-He solutions, should 
be used with high caution until a spectrum is secured. Among 
the other possible but unconfirmed white dwarfs, we derived the 
atmospheric parameters of LP 652-342 (probable non-member) 
in a similar fashion. For 1RXSJ062052.2+ 132436, we find that 
the photometry is consistent with a very hot object, but it is in a 
regime where r e jf can not be estimated with a reasonable accu- 
racy, and therefore it is omitted from Table 1 . 
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Fig. 2: Johnson BVRI and 2MASS JHK S photometry of 
WD 0120-024 (error bars) fitted by average model fluxes (filled 
points) for a pure-hydrogen atmosphere. The best fit parameters 
are given in the figure, and the corresponding predicted spectrum 
is also shown (solid line). 
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Fig. 3: Similar to Fig. |2] but for the BV and 2MASS JH fit of 
LP 246-14 using pure-hydrogen models (left panel) and pure- 
helium models (right panel). 



4. Spectroscopic and Kinematic distances 

The set of atmospheric parameters that we have compiled in 
Table 1 can be used, along with model atmospheres, to predict 
absolute magnitudes for any photometric band. By combining 
these values with observed magnitudes, we can compute the so- 
called spectroscopic distances (t/ spe ctiofl This is independent of 
the kinematic distances derived in Paper II to which they can be 
compared. Our photometric data set, including 2MASS JHKs, 
CMC14 r' and Johnson BV can be used to predict as much as 6 
alternative d spectro values. In Fig.|4]we present, for all observed 
magnitudes, the ratio of the spectroscopic and kinematic dis- 
tances for objects in Table 1. In the computation of the uncer- 



2 Note that the distances could also be determined with this method, 
but the technique described in Sect. 4 is more precise. 



3 Note that in a few cases, photometric instead of spectroscopic fits 
were used to determine the atmospheric parameters, but we use the term 
spectroscopic distance in this work for clarity. 
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Table 1 : Atmospheric parameters of Hyades white dwarf candidates 



Star No. 


WD name 


Other name 




SpT 


T e g 


log 9 


M/M Q 


M v 


T 


dist 


Ref. fl 












(K) 


(cm s~ 2 ) 




(mag) 


(yr) 


(pc) 




1 


0347+171 


V471 Tau 




DA+dK 


34100 (100) 


8.25 (0.05) 


0.84 (0.04) 


10.20 (0.12) 


7.07 (0.13) 


- 


1 


2 


0352+096 


HZ4 




DA 


14670 (380) 


8.30 (0.05) 


0.80 (0.03) 


11.76 (0.09) 


8.53 (0.05) 


36 


2 


3 


0406+169 


LB 227 




DA 


15810 (290) 


8.38 (0.05) 


0.85 (0.03) 


11.76 (0.09) 


8.50 (0.04) 


52 


2 


4 


0418+137 


HD 27483 




DA+dF 


21410(210) 


8.00 (0.50) 


0.63 


10.65 


7.70 


- 


3 


5 


0421+162 


LP 415-46 




DA 


20010 (320) 


8.13(0.05) 


0.70 (0.03) 


10.96 (0.08) 


7.97 (0.06) 


47 


2 


6 


0425+168 


LP 415-415 




DA 


25130 (380) 


8.12(0.05) 


0.71 (0.03) 


10.53 (0.08) 


7.49 (0.08) 


50 


2 


7 


0429+176 


HZ 9 




DA+dM 


17620 (350) 


8.02 (0.06) 


0.63 (0.04) 


11.02 (0.09) 


8.08 (0.06) 


39 


2A 


8 


0431+126 


HZ 7 




DA 


21890 (350) 


8.11 (0.05) 


0.69 (0.03) 


10.78 (0.08) 


7.78 (0.07) 


49 


2 


9 


0437+138 


LP 475-242 




DBA 


15120 (360) 


8.25 (0.09) 


0.74 (0.06) 


11.65 (0.15) 


8.47 (0.07) 


45 


4 


10 


0438+108 


HZ 14 




DA 


27540 (400) 


8.15 (0.05) 


0.73 (0.03) 


10.40 (0.08) 


7.30 (0.09) 


49 


2 


11 


0120-024 


LP 587-53 




DA 


5880 (80) 


8.12(0.16) 


0.66 (0.10) 


14.53 (0.24) 


9.51 (0.15) 


- 


5 


14 


0230+343 


GD 30 




DA 


15270 (250) 


7.98 (0.05) 


0.60 (0.03) 


11.21 (0.08) 


8.26 (0.04) 


91 


2 


15 


- 


LP 246-14 




- 


9680 (2600) 


8.00 (0.50) 


0.60 


12.31 


8.82 


- 


6 


16 


0259+378 


GD 38 




DA 


32040 (470) 


7.84 (0.05) 


0.57 (0.03) 


9.59 (0.09) 


6.89 (0.02) 


174 


2 


17 


0312+220 


GD 43 




DA 


18600 (300) 


7.89 (0.05) 


0.56 (0.03) 


10.74 (0.08) 


7.86 (0.06) 


105 


2 


18 


0339-035 


LP 653-26 




DA 


13000 (210) 


8.06 (0.05) 


0.65 (0.03) 


11.60(0.07) 


8.51 (0.04) 


52 


2 


19 


0348+339 


GD 52 




DA 


14820 (350) 


8.31 (0.05) 


0.80 (0.03) 


11.76 (0.09) 


8.52 (0.05) 


49 


2 


20 


- 


HS 0400+1451 




DA 


14620 (60) 


8.25 (0.01) 


0.76 (0.01) 


11.69 (0.02) 


8.50 (0.01) 


- 


7 


21 


0433+270 


LP 358-525 




DA 


5630 (100) 


8.06 (0.04) 


0.62 (0.02) 


14.65 (0.10) 


9.51 (0.06) 


18 


8 


22 


0437+122 


LP 475-249 




DA 


13260 (350) 


8.08 (0.08) 


0.66 (0.05) 


11.60 (0.12) 


8.50 (0.06) 


165 


2 


23 


0625+415 


GD74 




DA 


17610 (280) 


8.07 (0.05) 


0.66 (0.03) 


11.10(0.08) 


8.12(0.05) 


60 


2 


24 


0637+477 


GD 77 




DAP 


14650 (590) 


8.30 (0.06) 


0.80 (0.04) 


11.76 (0.12) 


8.53 (0.06) 


40 


2B 


25 


0641+438 


LP 205-27 




DA 


16400 (270) 


8.00 (0.05) 


0.62 (0.03) 


11.12 (0.08) 


8.17 (0.05) 


76 


2 


26 


0743+442 


GD 89 




DA 


15220 (350) 


8.47 (0.05) 


0.91 (0.03) 


11.98 (0.09) 


8.61 (0.05) 


39 


2 


27 


0816+376 


GD 90 




DAP 


11000 (630) 


8.00 (0.50) 


0.60 


11.87 


8.67 


- 


9 


28 


0233-083.1 


SDSS J023614.44- 


-080804.9 


DQA 


10000 (1000) 


8.00 (0.50) 


0.60 


12.18 


8.78 


- 


10 


29 


0300-083.1 


SDSS J030325.22- 


-080834.9 


DA 


11550 (200) 


8.55 (0.08) 


0.96 (0.05) 


12.62 (0.15) 


9.00 (0.09) 


- 


11 


30 


- 


LP 652-342 




- 


10940 (2840) 


8.00 (0.50) 


0.60 


11.88 


8.68 


- 


6 


31 


0533+322 


G98-18 




DA 


16970 (280) 


7.59 (0.05) 


0.43 (0.02) 


10.47 (0.08) 


8.09 (0.03) 


159 


2 


32 


0543+436 


G96-53 




DAZ 


8260 (130) 


7.46 (0.13) 


0.36 (0.05) 


12.19 (0.19) 


9.03 (0.07) 


99 


2 


33 


0557+237 


G104-10 




DA 


8360 (120) 


8.18(0.06) 


0.70 (0.04) 


13.15(0.11) 


9.10(0.05) 


56 


2 


35 


0758+208 


SDSS J0801 14.02+204340.3 


DA 


7640 (70) 


8.68 (0.20) 


1.03 (0.12) 


14.35 (0.38) 


9.56 (0.04) 




12 


36 


0816+387 


Glll-71 




DA+dM 


7700(110) 


8.07 (0.07) 


0.64 (0.04) 


13.31 (0.12) 


9.12 (0.05) 


45 


2A 


37 


0820+250.1 


SDSS J082346. 14+245345.6 


DA 


34050 (100) 


7.74 (0.02) 


0.53 (0.01) 


9.29 (0.03) 


6.80 (0.01) 




11 


38 


0231-054 


PHL 1358 




DA 


17470 (270) 


8.58 (0.05) 


0.98 (0.03) 


11.94 (0.09) 


8.52 (0.04) 


29 


2 


39 


0548-001 


G99-37 




DQP 


6130 (100) 


8.18 (0.04) 


0.69 (0.03) 


14.42 (0.10) 


9.56 (0.05) 


11 


8 


40 


0648+368 


GD 78 




DC 


5700(1000) 


8.00 (0.50) 


0.59 


14.51 


9.42 




13 



Notes. Star Nos. and categories are the same as in Paper II. Nos. 1 to 10 are the so-called classical Hyades white dwarfs. Nos. 1 1 to 27 are new candidates (including star Nos. 23 to 27 that are more 
than 20 pc away from the center in the Z direction). Nos. 28 to 37 are rated probable non-Hyades members. Nos. 38 to 40 are new candidates identified in Sect. 2 with a residual velocity upper limit 
of 8 km s -1 . 

References. (1) ISion et al] d2012h: (2) iGianninas et all J201 ll); (2A) same as previous but with a M-dwarf component removed; (2 B) same as previous b ut magnetic white dwarf fitted with 
non-magnetic model s; (3 ) iBurleigh et al.l dl998l); ( 4) iBergeron et al] d2011l); (5) photometric fit fo r r eff , this work, and mass from ISilvestri et alj d2001l ) ; (6) p hoto metric fit, this work; (7 ) 
iKoester et all d2009bh; (8) lGiammichele et alj d2012l) ; (9) lLimoges & Bergeron! d2010h ; lJordad dl993l) ; (lO) IVanlandingham et al] d2005l) ; (1 D lTremblav et"alTd201 lal) ; (12) lEisenstein et al] ( 12000) ; 
(13) lAngeletal.ldl981h . 
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Fig. 4: Ratio of the spectroscopic to kinematic distance for the Hyades white dwarf candidates identified by their star number from 
Table 1 (from No. 1 to 20 in upper panel and from No. 21 to 40 in the lower panel). Up to five <i sp ectro values were evaluated for each 
candidate using 2MASS JHK$ (red), CMC r' (cyan) and Johnson BV (blue) photometry (open points). Many B and V magnitudes 
have no published uncertainties. Also shown is the the mean weighted spectroscopic distance from 2MASS JH, CMC r' and V 
filters (black filled points). The data points are shifted horizontally for clarity. Objects with a cross above the star No. are binaries 
(Nos. 1, 4, 7, 12 and 13) or have no atmospheric parameters (No. 34). For several objects, the ratio is actually higher than the scale 
of the plot, and is represented by a lower limit symbol. The dotted horizontal line represents a perfect match between both distance 
estimates. 



tainties, we account for the contribution of the d^ n prediction 
(depending only on the observed proper motion and sky posi- 
tion), the atmospheric parameters and the observed magnitude. 
We note that unresolved binaries are omitted from Fig. [4] since 
the observed magnitudes are all either mildly or strongly con- 
taminated by the companion. 

It is clear that the different d spectIO determinations for a same 
object are not completely independent, since they are based 
on the same set of atmospheric parameters. Nevertheless, by 
looking at the results between different filters, one can get an 
overview of the quality of the photometric data sets. It is seen 
that the scatter is in general low between the different filters, 
implying that our photometric data sets are fairly homogeneous 
and well calibrated. We nevertheless note that one Kg magni- 
tude (No. 16) and one H magnitud e (No. 20) in Fig. Blare signif- 
icantly different to other values. iTremblav & Bergeronl |2007) 
have shown that one must be careful with spurious magnitude 
values in 2MASS close to the survey detection limit, and these 
are unlikely to be the signature of real near-infrared spectral fea- 
tures. Furthermore, B and V magnitudes are drawn from differ- 
ent sources with inhomogeneous accuracy, and often even have 
no quoted uncertainties Based on these facts, we decided to 



4 Our results indicate that B and V magnitudes are nevertheless 

in fairly good agreement with the other data sets, and that the of- 

ten used typical average uncertainty of 0.05 mag in V (see, e.g., 
ITremblav & Bergeronll2008h appears adequate. 



compute a weighted mean d spectl0 value using only 2MASS JH, 
CMC r' and V magnitudes (when the uncertainties are givenfl 
In Fig. H] we show the comparison between the mean d spectIO 
and d\d„ determinations for each candidate. These quantities are 
also given in Table 2, along with the number of filters involved 
in the d spectIO mean. Note that at this point, we leave out from 
our calculations and discussion the objects labelled as probable 
non-members in Paper II (Nos. 28-37), since they are clearly all 
background objects according to Fig. |4] 

In the next section we will discuss about the Hyades mem- 
bership of new candidates based on the results of Fig. [4] 
Beforehand, there is some interesting information that can be 
extracted from the classical members. The non-binary Hyades 
white dwarfs (Nos. 2-3, 5-6 and 8-10) can serve as a bench- 
mark to evaluate the precision of our analysis. The spread due 
to the photometric uncertainties has already been evaluated from 
Fig- 111 It is now seen from the same figure that the mean d sp&cao 
values derived in this work are in good agreement, roughly 
within the uncertainties, with d^ n values. The spread around the 
expected unity ratio can be used to learn about the combined 
uncertainties due to spectroscopic parameter determinations and 
<4in predictions. The good agreement suggests that the uncertain- 
ties in this work are relatively well estimated and that our method 
should be fairly robust to argue on the status of the Hyades mem- 

5 For stars No. 22 and 40, we assume a V uncertainty of 0.2 since no 
2MASS or CMC measurements are available. 
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bership of new candidates. Furthermore, we note that the sin- 
gle classical members have a mean mass of 0.75 M©, signifi- 
cantly higher than the value of 0.638 M , with a dispersion o f 
0.145 M e , found for field white dwarfs dGianninas et al.ll201 lh . 
Also, their cooling ages are in the range 10-340 Myr, which 
is consistent with the young age of the Hyades, though signif- 
icantly smaller than the estimated a ge of the cluster of 625 Myr. 
According to iDobbie et alj d2006l) . this corresponds to main- 
sequence progenitors in the range of ~2.5-4 M . Therefore, 
more massive and slightly older white dwarfs, hence fainter and 
cooler objects, could be expected among the new candidates. 

One uncertainty that is often difficult to evaluate, and not in- 
cluded in Fig. |4] is the precision of the model atmospheres. The 
classical Hyades white dwarfs, however, can help us for the eval- 
uation of this uncertainty. Since any possible inaccuracies in the 
model atmospheres are expected to cause a systematic offset of 
the atmospheric parameters, they would also cause systematic 
shifts of the points in Fig|4] Other uncertainties discussed so far 
are much more likely to cause a spread around a mean value. Our 
results show that there is no need to account for an additional off- 
set, and therefore we conclude that the atmospheric parameters 
are sufficiently accurate on the average. This is a fairly precise 



validity verification of the physics in the Tremblav & Bergeron 
(2009) DA models, and to a lesser degree the Ber geron et al 



(2011) DB models (since only star No. 9 is helium-rich), which 
were adopted to predict the atmospheric parameters of classical 
members. Reversing this statement, it illustrates that the method 
of using model atmospheres and spectroscopic observations to 
predict t/ S pectro values, which has become a standard method in 
recent years (see, e.g ., iHolberg et all 120081; IDobbie et al1l2009t 
Willi ams et al.ll200"9l) . is fairly well calibrated and precise on av- 
erage. One still has to be cautious about the precision of ob- 
served magnitudes used in this process. 



5. Membership status of new Hyades candidates 

We have tested our t/ spe ctro determinations with the classical 
Hyades, and we expect these values to be equally precise for 
the new aspirants. In this section, we review the status of the 17 
new candidates (Nos. 1 1-27) identified in Paper II and the 3 new 
candidates (Nos. 38-40) identified in this work. For that purpose, 
we have computed in Table 2 the new distances between the ob- 
jects and the Hyades center (c/Hyades center) based on the af S pectro 
values. It is clear that all but one (No. 20) of the new candidates 
are outside of the tidal radius. 

Another quantity in Table 2 derived from the t/ spe ctio calcu- 
lation is the residual tangential velocity |y ta ngl with respect to 
the expected motion of the star towards the cluster convergent 
point. In Fig. [5] the residual velocity of the candidates is plotted 
versus the distance to the cluster center. Our updated |u t angl val- 
ues are by definition higher than those found in Paper II, where 
the difference in the proper motion of candidates and that of the 
Hyades where minimized on a plane perpendicular to the line-of- 
sight, in order to derive the kinematic distances. Also included 
in Fig. [5] are the two unresolved binaries for which we have no 
atmospheric parameters (Nos. 12 and 13), and we use their kine- 
matic parameters in order to have the complete picture of possi- 
ble members. Because the companion of these candidates follow 
very well the Hyades main sequence in the absolute magnitude 
vs. color diagram for near-infrared colors (see Fig. 4 of Paper I), 
their kinematic distances are expected to be fairly accurate. 

According to Fig. |U the c/ S pecti-o and dkm determinations do 
not agree within lcr for most candidates, which would imply 
that their membership is doubtful. On the other hand, if we al- 
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Fig. 5: Residual tangential velocity \v tang \ with respect to the 
Hyades group motion as a function of the distance to the clus- 
ter center, for classical members (blue points and error bars) and 
new candidates with an age lower (black points and error bars) 
and higher (red points and error bars) than the Hyades age. Also 
shown are the two candidates in binaries (Nos. 12 and 13, cyan 
points in the figure) for which we have no spectroscopic param- 
eters. Instead we rely on the kinematic values, providing lower 
limits on the residual velocities. The error bars are not rms er- 
rors but refer to the minimum and maximum values where the 
distance of a given white dwarf is varying within d spectro + lcr. 
The horizontal dotted line represents the 8 km s _I residual ve- 
locity upper limit for the candidate selection. 



low them to have residual velocities of a few kilometers per sec- 
ond, they could still be valid unbound constituents. As a conse- 
quence, the residual velocities of Fig. [5] are better suited to look 
at the membership status of new candidates. The residual veloc- 
ities are of the order of ~1 km s for the classical members, and 
one could in principle use this characteristic value for a mem- 
bership confirmation. However, Paper I has shown that unbound 
main-sequence stars had a significantly higher residual velocity 
dispersion than the bound members, and we can expect a similar 
behaviour for white dwarfs. Then, a higher residual velocity can 
be allowed to judge the status of the new candidates. 

It is already clear that 6 of the candidates (star Nos. 14, 16, 
17, 18, 22 and 25) must be ruled out as Hyades components, 
since their updated distances imply residual velocities of more 
than 20 km s _1 . Indeed, it is obvious that a cluster candidate can 
not have a residual velocity of the order of the velocity disper- 
sion of fi eld white dwarfs, which is ~30 km s _1 in each spatial 
direction (ISion et al.ll2009l) . 

Before we look at the remaining candidates, it is desirable to 
understand why the above 6 background objects contaminated 
the sample of Paper II. From the kinematic selection alone, it 
is clear that background, high velocity objects, could pollute 
the sample. However, these background stars should have been 
removed in the photometric selection, where more distant ob- 
jects appear subluminous compared to true Hyades members 
(see Fig. [T]i. In Fig. [6] we compare the observed and predicted 
B - V colors, the latter computed from model spectra with the 
atmospheric parameters from Table 1. We conclude that only 
two objects (Nos. 14 and 17) have been misidentified as can- 
didates since their published B — V colors erroneously put them 
on the white dwarf sequences. For the other four candidates re- 
jected as members, we conclude that they just happened to be 
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Table 2: Distances of Hyades white dwarf candidates 



Star No. 


Name 


4dn 




Filters" 


i lb 

If tang 1 


/ * 

"Hyades center 


Status 1, 






(pc) 


(pc) 




(km s ') 


(pc) 




2 


WD 0352+096 


35.1 (0.7) 


35.4 (1.4) 


4 


1.0 (0.9, 1.7) 


13.3 (12.2, 14.4) 


classical 


3 


WD 0406+169 


51.8 (1.6) 


51.5 (2.0) 


4 


0.5 (0.4, 1.2) 


6.5 ( 5.0, 8.2) 


classical 


5 


WD 0421+162 


44.1 (1.4) 


46.0 (1.7) 


4 


1.4 (0.9, 2.2) 


1.0 ( 0.9, 2.2) 


classical 


6 


WD 0425+168 


47.8 (2.1) 


49.9 (2.0) 


3 


1.2 (0.6, 2.2) 


3.6 ( 1.6, 5.6) 


classical 


8 


WD 0431 + 126 


45.7 (0.4) 


48.6 (1.9) 


4 


1.5 (0.6, 2.4) 


4.2 ( 3.4, 5.5) 


classical 


9 


WD 0437+138 


46.9 (1.8) 


45.4 (3.2) 


4 


1.0 (0.7, 2.2) 


3.3 ( 3.2, 5.2) 


classical 


10 


WD 0438+108 


46.1 (1.1) 


48.6 (1.9) 


4 


1.2 (0.4, 2.0) 


5.8 ( 5.3, 6.9) 


classical 


11 


WD 0120-024 


39.1 (0.8) 


42.0 (5.0) 


2 


3.5 (0.8, 9.1) 


37.1 (35.9, 39.0) 


old WD 


14 


WD 0230+343 


37.9 (0.6) 


91.1 (4.5) 


4 


56.0 (51.3, 60.8) 


56.8 (52.8, 60.9) 


non-member 


15 


LP 246-14 


35.1 (0.5) 


69.1 (33.4) 


4 


37.9 (2.5, 75.0) 


36.8 (23.3, 66.3) 


old WD 


16 


WD 0259+378 


66.2 (1.8) 


172.3 (6.9) 


3 


61.9 (57.9, 65.8) 


133.4 (126.6, 140.2) 


non-member 


17 


WD 0312+220 


44.4 (1.0) 


101.0 (5.6) 


4 


45.3 (40.9, 49.7) 


58.7 (53.3, 64.1) 


non-member 


18 


WD 0339-035 


28.5 (0.4) 


53.5 (2.1) 


3 


27.5 (25.4, 29.8) 


21.2 (20.2, 22.4) 


non-member 


19 


WD 0348+339 


38.5 (0.7) 


47.1 (2.1) 


3 


7.3 (5.6, 9.1) 


15.6 (15.3, 16.2) 


candidate 


20 


HS 0400+1451 


38.5 (2.5) 


45.2 (1.3) 


3 


5.1 (4.2, 6.0) 


5.1 ( 5.0, 5.4) 


candidate 


21 


WD 0433+270 


20.1 (0.2) 


17.3 (0.7) 


4 


3.6 (2.7, 4.5) 


29.5 (28.8, 30.2) 


old WD 


22 


WD 0437+122 


67.1 (7.7) 


166.1 (17.9) 


1 


31.7 (26.0, 37.4) 


120.0 (102.2, 137.9) 


non-member 


23 


WD 0625+415 


47.6 (1.4) 


58.0 (2.4) 


3 


6.6 (5.5, 7.9) 


34.0 (32.6, 35.5) 


candidate 


24 


WD 0637+477 


36.2 (0.7) 


39.3 (2.1) 


3 


3.0(1.7, 4.6) 


30.9 (30.6, 31.2) 


candidate 


25 


WD 0641+438 


41.8 (1.1) 


71.2 (3.1) 


3 


19.9(17.9,21.9) 


46.2 (43.8, 48.6) 


non-member 


26 


WD 0743+442 


33.6 (0.6) 


38.1 (1.7) 


4 


4.3 (2.9, 5.8) 


36.4 (36.0, 36.8) 


candidate 


27 


WD 0816+376 


38.8 (0.9) 


54.5 (18.8) 


3 


11.8(0.9, 25.9) 


46.8 (38.8, 59.9) 


uncertain 


38 


WD 0231-054 


32.8 (0.4) 


27.2(1.3) 


3 


8.6 (7.6, 9.7) 


29.0 (28.6, 29.5) 


candidate 


39 


WD 0548-001 


9.4 (0.2) 


11.1 (0.4) 


3 


6.8(6.4,7.1) 


36.7 (36.4, 37.1) 


old WD 


40 


WD 0648+368 


36.2(1.0) 


31.8 (17.4) 


1 


5.7 (5.6, 14.1) 


28.6 (28.2, 36.0) 


old WD 



(a) Number of photometric measurements used to compute the mean spectroscopic distance. 

{b) Numbers in parentheses are not rms errors but refer to the minimum and maximum values where the distance of a given white dwarf is varying 
within rf sp ec(ro ± lcr 

(c) Summary of the status of the Hyades membership as discussed in this paper, including the classical members and the new candidates that 
are still valid candidates or that were rejected as non-members following our analysis. We also identify old white dwarfs for which the age is 
incompatible with a Hyades origin. Finally, one object has uncertain atmospheric parameters and we can not rule it out as a member. Note that 
binaries (three classical members, Nos. 1, 4 and 7, and two new candidates, Nos. 12 and 13) are excluded from this table. 

difficult, from an absolute magnitude vs. color diagram alone, 
to conclude if a kinematic candidate is a Hyades member, or a 
hotter and more distant background white dwarf. The method of 
comparing kinematic and spectroscopic predictions used in this 
work is then necessary to obtain a clean sample. One might be 
intrigued by the fact that we find 6 background objects (up to 172 
pc from the Sun), and no foreground objects. The main reason 
seems that the volume surveyed in this work in front of the clus- 
ter is only -1/50 of the volume beyond the cluster up to the most 
distant background object. This volume effect may not be im- 
portant since catalogues of white dwarfs are far from complete 
at distances greater than the Hyades center. On the other hand, 
all of our background objects are hot DAs, and these objec ts are 
fairly rare up to a radius of 20 pc ( Giammic hele et al1l2012l) and 
are more probably found at larger distances. 

There remain 14 objects, including 2 binaries with no spec- 
troscopic parameters, that are still possible members. Five can- 
didates (Nos. 11, 15, 21, 39 and 40), however, are significantly 
older than the age of the cluster, including one object (No. 15) 
that is not spectroscopically observed and hence confirmed as 
a white dwarf. Furthermore, we have one magnetic object with 
imprecise atmospheric parameters (No. 27). According to Fig. [5] 
most of the remaining candidates have significant residual veloc- 
ities of the order of 5 km s _1 and lie far away outside of the tidal 
radius. These candidates, if they are of Hyades origin, are not 
clearly linked to the classical members in terms of their clus- 
ter center distances and residual velocities, but could be escaped 
former members. 



> 
i 



> 
I 

a 




1 2 3 4 5 6 7 



10 11 12 13 14 15 16 1? 18 19 20 21 22 23 24 25 26 27 



Star number 



Fig. 6: Difference in the observed and predicted B-V colors for 
the Hyades white dwarf candidates identified by their star num- 
ber. The dotted horizontal line represents a perfect match be- 
tween both quantities. Objects with a cross above the star No. 
are binaries. 



within the theoretical sequences by chance, but with a vertical 
offset. Indeed, for hot white dwarfs, the absolute V magnitude 
is rather insensitive to the B-V color as can been seen from 
Fig. Q] More precisely, for a given color, a range of My of as 
much as two magnitudes is possible. It implies that it is very 
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In principle, the diagram of Fig. [5] should be fairly complete 
up to 8 km s~ 1 , although it must be kept in mind that not all white 
dwarfs around the Hyades have been identified in white dwarf 
catalogues. The goals of the next sections will be to understand 
the contamination from field white dwarfs and the completeness 
of the white dwarf catalogues. 

5.1. Field white dwarfs contamination 

We made Monte Carlo simulations to estimate the number of 
field white dwarfs that could lie within a radius of 40 pc of the 
cluster and mimic, by chance, the Hyades motion, hence cre- 
ating false positive detections. As input, we need the luminosity 
function of white dwarfs within ~100 pc of the Sun. Such a func- 
tion is not known, but, b ased on improved atmosp heric parame- 
ters of degenerate stars. iGiammichele et al.l C2012) published the 
luminosity function of the relatively complete sample of white 
dwarfs within a 20 pc radius from the Sun. In principle, it would 
be appropriate to extrapolate their results to the Hyades region. 
However, the hot end of the luminosity function, which is of sig- 
nificant interest for the young cluster, suffers from low number 
statistics in the local sample of white d warfs. Conseque ntly, we 
rely on the SDSS luminosity function (lHarris et al.l2006l) . which 
includes a much larger number of hot objects, but is incomplete 
by as much as 50% and complex completen ess corrections must 
be ac counted for. Nevertheless, Fig. 22 of IGiammichele et alj 
d2012l) shows that the local sample and SDSS luminosity func- 
tions agree relatively well, and in particular, the derived total 
space densities are very similar. Therefore, we simulated white 
dwarf populations weighted by the SDSS luminosity function 
in a sphere of 40 pc radius centered on the Hyades. The objects 
are randomly distributed in the galactic X, Y, Z coordinates, with 
axes pointing to the Galactic Center (X), the direction of galactic 
rotation (Y), and the North Galactic Pole (Z). The U, V, W ve- 
locities are distributed according to a Gaussia n weight function 
using the observed velocity distribution from lSion et al.l ([2009) 
for the local 20 pc sample. They find mean velocities, with re- 
spect to the Sun, of [<U>, <V>, <W>] = [-0.9, -21.5, -5.2] 
km s and dispersions of [o~u, cry, <Tw] = [35.3, 33.1, 28.1] 
km s" 1 . To derive the mean velocity and the velocity dispersion, 
ISion et al ] (12009b assumed zero radial velocity, which, at least, 
sets a lower limit to the velocity dispersion. 

Our 100 simulations yield on average 1233 white dwarfs 
within 40 pc of the cluster center. Nonetheless, we do not ex- 
pect that many of these objects will have velocities close to the 
Hyades group motion of [-41.1, -19.2, -1.4] km s _1 (Paper I). 
On average, we find that 2.6 objects have 3D velocities that are 
consistent with the Hyades motion within 8 km s . However, 
3D velocities can not be measured in most cases due to the lack 
of radial velocities. Hence we also measured the residual 2D ve- 
locities of our simulated objects with respect to the Hyades mo- 
tion perpendicular to the line-of-sight. In that case we find on av- 
erage 29 objects that mimic the cluster tangential motion within 
8 km s . We show one of our simulations in the upper panel of 
Fig. [71 with the relative 2D velocity as a function of the cluster 
center distance. The lack of objects close to the center is simply 
a volume effect. We note that cool white dwarfs dominate the 
luminosity function and also the local white dwarf sample, and 
consequently 88% of our simulated candidates are cool objects 
(My > 12, red points on the figure), for which the age, above 
~500 Myr, would be inconsistent with a Hyades membership. In 
other words, only about 3 of the contaminants would have an age 
consistent with the Hyades age. The summary of the predicted 
number of polluting objects is given in Table 3. Finally, it should 



be realized that due to the uncertainties in the observed luminos- 
ity functions and velocity distributions, the numbers given in this 
section are only accurate to about the 10% level. 
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Fig. 7: Top: Residual tangential velocity with respect to the 
Hyades group motion, as a function of the distance to the clus- 
ter center, for a selected Monte Carlo simulation of field white 
dwarfs as described in the text. The red points are cool white 
dwarfs with My > 12. The horizontal dotted line represents the 
residual velocity upper limit of 8 km s _! used in Section 5.1. 
Bottom: Same simulation as the top panel, but corrected for the 
completeness of the WDC (see Section 5.3). 



The selection criteria of candidates in Paper II include a 
cross-match with the White Dwarf Calatog (WDC) as well as 
with a few other sources. These catalogues suffer from incom- 
pleteness and we have to take this into account in our Monte 
Carlo simulations. In the following, we therefore investigate the 
incompleteness of the WDC to determine a realistic estimate of 
the field star contamination in our sample. 



Table 3: Contamination from field white dwarfs 



Residual velocity (< 8 km s ') 


Total 


M v < 12 


M v > 12 


uncorrected for completeness 


3D 


2.6 


0.30 


2.3 


2D 


29 


3.3 


26 


corrected for completeness 


3D 


0.34 


0.21 


0.13 


2D 


3.9 


2.4 


1.5 



Notes. Given a maximum cluster center distance of 40 pc. 
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5.2. Completeness corrections 

The V magnitude distribution of the simulated white dwarfs 
from the previous section is shown in Fig. [8] (top panel, bold 
histogram). Since the local sample of white dwarfs is only com- 
plete up to about 20 pc and that the objects in our simulations 
are on average twice as distant, it is already obvious that many 
of them may be missing in white dwarf catalogues. One impor- 
tant thing to notice, however, is that all of the faint objects with 
V > 17 in our simulations are cool white dwarfs with My > 12. 
Indeed, for white dwarfs with My < 12, i.e. those with an age 
roughly consistent with the age of the Hyades, the bottom panel 
of Fig.[8]shows that they are always relatively bright (V < 17). 




10 11 12 13 14 15 16 17 18 19 20 21 22 



V (mag) 



Fig. 8: Top: Distribution of V magnitudes for the simulations of 
field white dwarfs discussed in Section 5.2 (bold histogram). 
The distribution has then been folded with our estimated com- 
pleteness fraction of the WDC (red hatched histogram). Bottom: 
Same as top panel but only for objects with My < 12. 



is fitted well by the expected relation. We computed the com- 
pleteness at V > 14.5 from the ratio of observed to predicted 
objects. Afterwards, we folded this completeness function with 
the distribution of V magnitudes in our field white dwarfs simu- 
lations presented in Fig. 10 (red dashed histograms). We derive 
from the new total number of objects that only 14% of the field 
white dwarfs within a radius of 40 pc to the Hyades center are 
found in the WDC. This fraction increases to 72% if we only 
account for young white dwarfs with My < 12, which implies 
that the completeness for old stars above this threshold is only 
about 6%. To find the number of field white dwarfs mimicking 
Hyades members, we simply have to multiply the completeness 
fractions to the numbers of the previous section. The result is 
given in Table 3, and, as an example, in Fig. [7] (bottom panel) 
where we modified the simulation by allowing objects to ran- 
domly disappear in accordance with the completeness correction 
as a function of magnitude. Hence, we predict that only 0.3 field 
white dwarfs, on average, will be found within 8 km s _I from 
the 3D cluster velocity. Furthermore, 3.9 objects, among those 
2.4 with My < 12 will mimic the 2D velocity within the same 
range. 



tso 
o 




19 20 



Fig. 9: Logarithm of the number of stars in the WDC as a func- 
tion of their observed visual magnitude, in bins of 0.5 mag. The 
error bars are from number statistics. A linear relation assuming 
a volume complete sample (dotted red line) as been fitted to the 
data in the range 12.5 < V < 14.5. 



To estimate the completeness of the White Dwarf Catalog, 
we used the following procedure. We computed the V mag- 
nitude histogram for all objects in the catalog which is pre- 
sented in Fig. [9] in terms of the logarithm of the number. We 
neglected all objects first discovered in the SDSS. Removing 
these objects seems justified to estimate the contamination of 
the Hyades sample, since all candidates in Table 2 were known 
before the SDSS era, and SDSS only partly covers the Hyades 
area. For uniform stellar density, a linear relation, with a slope 
of 0.6, is expected if a sample is volume complete. We can see 
that the observed distribution in the range 12.5 < V < 14.5 



6 We have neglected magnitudes taken prior to 1972, photographic 
magnitudes, and those flagged as uncertain. 



Although the relative larger number of white dwarfs at V < 
12.5 in Fig. [9] can be explained by statistical fluctuations, they 
can also indicate that our assumption of a completeness of 100% 
up to V — 14.5 is rather optimisticQ A larger completeness cor- 
rection provides a lower number of field white dwarfs contami- 
nating our sample of Hyades candidates. Therefore the contam- 
ination estimates above can be considered as conservative. 



7 Note that for V < 11, most objects are white d warfs in Siriu s-like 
systems, for which the completeness is rather low (Holbera 2009). 
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6. Discussion 

Our analysis has so far reviewed the status of the ten classical 
Hyades white dwarfs and provided improved distances based on 
the most recent determinations of the atmospheric parameters. 
Then, we updated the residual velocities with respect to the bulk 
motion of the cluster. We have furthermore studied 20 new can- 
didates identified in Paper II and in this work. We found that all 
but one were outside of the tidal radius, and six of them were 
clearly background objects. There remain 14 objects for which 
the Hyades membership status is still possible. A comparison of 
the distribution of the observed candidates versus the simulated 
field white dwarfs in Fig. [5] and Fig. [7] (bottom panel, with the 
completeness correction) can serve to conclude about the global 
status of those remaining candidates. 

We will first draw our attention on objects that have an age 
lower than the cluster age (black points in Fig. [5}, which are 
the most logical candidates. There are six such white dwarfs, 
although this excludes the two binary candidates for which the 
cooling age is not known (cyan points). Also, the magnetic white 
dwarf WD 0816+376 is at 47 pc from the cluster center accord- 
ing to Table 2, but given the imprecise atmospheric parameters, 
we can not exclude that it is actually closer (see Table 2). In com- 
parison, the number of contaminants in this age group is 2.4 or 
3.3, including or not the completeness correction, respectively. 
We find that the probability that six objects are field stars is 4% 
or 12% for a Poisson distribution including or not the complete- 
ness correction, respectively. Compared to the field, there is an 
overdensity of white dwarfs of ages younger than the Hyades. 
We must remind, however, the difficulty of estimating the uncer- 
tainty on the number of contaminants, although we believe that 
our proposed number is conservative. 

If we look closer at the properties of the 6 most promising 
candidates (identified in Table 2 by their candidate status), we 
can see from Table 1 that their masses are in the range of 0.66 
< M Q < 0.98, with a mean of 0.82 M Q , i.e. much higher than the 
mean mass of field white dwarfs. The probability of 6 field stars 
having such a high mean mass, assuming the v olume-complete 
mass distribution of Giammi chele et all d2012l) . is only ~0.1%. 
Once again, these 6 objects stand out from the field white dwarfs. 

Secondly, we can look at older stars (My > 12, red points 
in the Fig. [5] and [7]) for which we have four candidates. For the 
time being, No. 15 remains an additional candidate, given its un- 
certain nature and imprecise atmospheric parameters. From our 
simulations, we expect between 1.5 and 26 old objects to con- 
taminate the sample, considering or not the completeness correc- 
tion, respectively. Since the completeness correction is so high, 
it is hard to be confident about the expected number of contami- 
nants. Taking literally the mean of 1.5 contaminants, the proba- 
bility to find 4 is only 7%. However, one would expect that these 
are the most obvious contaminants, since their ages are inconsis- 
tent with the Hyades age. 

Nevertheless, an important problem arises if one concludes 
too fast that the old white dwarfs are contaminants. Indeed, two 
of these candidates have radial velocities t hat are in agreement 
with the Hyades m otion (WD 0120-024; ISilvestri et alJl200lL 
and WD 0433+270; lZuckerman et aU 2003). which would make 
them strong candidates if it was not for their age of about 4 
Gyr (see Table 1). These two objects are in common proper 
motion pairs with main-sequen ce stars consid e red as of Hyades 
origin. We note, however, that ISilvestri et a il (1200 lh also mea- 
sured a radial velocity for WD 0433+270, which is incompat- 
ible with the value found by IZuckerman et alj ((2003), and also 
incompatible with the radial velocity of its companion, since the 



mass, hence the gravitational redshift correction, is known for 
this white dwarf. Our simulations show that we expect to find, 
on average, 2.3 old objects that mimic the 3D Hyades velocities 
within 8 km s _1 , although when accounting for the completeness 
correction, the number drops to about 0. 1 . The fact that we find 
two of them may be statistically improbable, but not impossible 
considering low number statistics. If we instead infer that they 
are not simply field stars that accidentally have the 3D Hyades 
motion, different scenarios can be conceived. Perhaps the most 
promising scenario is that these objects could be members of 
the Hyades supercluster cold kinematic stream which is known 
to have components that are significantl y older than the Hyades 
(IChereul et alJll999l : lFamaev et al.ll2007l) . Other possibilities are 
that they could have been captured by th e Hyades cluster, or they 
could be younger Fe-core white dwarfs dCatalan et al.ll2008h . al- 
though in the latter case, the authors claim that the white dwarfs 
would still be too old (~1 Gyr). 

As a consequence of the discussion above, further obser- 
vations are needed to constrain the membership of the Hyades 
candidates. Accurate trigonometric parallaxes, for instance from 
Gaia, would provide distances that are independent of d spectm 
and therefore can eliminate possible uncertainties due to the at- 
mospheric parameters. This alone, however, is not sufficient to 
conclude about their membership, since their true 3D veloci- 
ties are not known. Consequently, accurate radial velocity mea- 
surements are indispensable to confirm or rule out which of the 
candidates is of Hyades origin. If consistent with the predicted 
Hyades radial velocity, the probability of contamination is very 
low (see Table 3). Due to the extreme broadening of the spectral 
lines, the spectroscopic measurement of the radial velocity of 
white dwarfs is extremely difficult to obtain, a nd only possible 
with the help of the narrow NLTE core s of H a dGreen stein et al. 
U977l:lKoesterll987HFalcon et alJ2010h . Since this measurement 
must be corrected for the gravitational redshift, typically ~70 
km s _1 , such a determination is only possible if we have a good 
knowledge of the mass. For instance, star No. 2 has a radial ve- 
locity measurement that is not in agreement with the Hyades mo- 
tion (see Paper II), even though it is still considered as a classi- 
cal member. It is more likely, in this case, that the uncertainty of 
the radial velocity measurement is underestimated. One should 
therefore be cautious about false negative conclusions from ra- 
dial velocities alone. 

In principle, there is a completely different way to mea- 
sure the radial velocity astrometrically by determining the 
perspective a cceleration due to the changing dista nce of the 
white dwarfs (IOortlll932tlGatewood & Russeliri974 with Gaia 
dAngla da-Escu de & Debesll2010h . In this approach no correction 
for gravitational redshift is needed. Assuming that a 3cr devia- 
tion of the proper motion during the 5 years mission is required 
and taking into account that th e Gaia proper-mo tion standard er- 
ror depends on the magnitude dde Bruiinell20Tl . one concludes 
that only bright (V < 13) white dwarfs closer than 5 pc allow 
an estimate of the radial velocity with the needed precision (de 
Bruijne, priv. comm.). Hence we can not rely on this method for 
the Hyades. 

6.1. Implications in terms of the Hyades evolution 

This paper has performed a spectroscopic check of the 27 
white dwarfs found in Paper II as possible Hyades members. 
Compared to Paper II, we tried to reveal more kinematic candi- 
dates by softening the velocity criterion. Nevertheless, no more 
than 10 white dwarfs have been found within the present-day 
tidal radius (9 pc) of the cluster, including one new candidate 
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(No. 20). Also, no objects with an age between 340 and 600 Myr 
are recovered in this radius range. Here we claim that no more 
single white dwarfs of ages lower than the Hyades age will be 
detected within the tidal radius. Outside the tidal radius and up 
to 40 pc, between 6 and 9 white dwarfs with a possible Hyades 
origin are found. Among those, about three are expected to be 
contaminants, although a few hot degenerate candidates could 
also be missing due to the incompleteness of the white dwarf 
catalogues. Even among the 6 most promising new candidates, 
all have an age lower than 340 Myr, except for No. 26, with an 
age of 41 Myr. In summary, we fa ce again the long-standing 
problem (Tins levlll974t IWeidemannl 1 1 9771) that not enough old 
white dwarfs are associated with the cluster in comparison to 
predictions, in particular co mpared to the mos t recent N-body 
Hyades cluster simulations of Ernst et al] d201 ll) . 

The simulations of Erns t et all (1201 ll) use the present-day 
mass distribution of stars within 9 pc of the cluster center from 
Paper I among inputs (see their Table 1) to simulate the for- 
mation and evolution of the cluster. Their best fit simulations 
ultimately predict the number, mass, and age of white dwarf 
members. They forecast that 73 remnants were formed in the 
cluster so far, including five black holes, 13 neutron stars and 
55 white dwarfs. However, many of these remnants have long 
ago evaporated from the cluster. In the simulations, the newly 
formed white dwarfs get a kick drawn from a Maxwellian distri- 
bution with a ID velocity dispersion of 5 km s _1 , corresponding 
to a 3D kick with a mean of ~8 km s . The present-day con- 
sequence of these kicks is that 40% (22 objects) of the white 
dwarfs remain within the present-day tidal radius and are dy- 
namically cooled down to show a velocity dispersion of only 
~0.5 km s 1 . It is worth noting that only ~ 10 among these bound 
white dwarfs have ages smaller than 340 Myr (Ernst et al., pri- 
vate communication). On the other hand, the kicks have removed 
33 white dwarfs from the cluster. At present, five lie in the range 
9 < ^center < 30 pc, and others are more than 30 pc away from 
the center. The authors are cautious to warn that the final number 
of tidally bound white dwarfs depends critically on the kick ve- 
locity. However, there are actually few observational constraints 
to confirm if white dwarfs obtain such kicks from asymmetri- 
cal mass l oss in the planetary nebula ejection, an d how strong it 
would be dDavis et a l. 2008; Freg eau et alj |2009). We also note 
that possible encounters with giant molecular clouds or spiral 
arms are neglected in the simulations. 



Compared to the N-body simulations of lErnst et alj (120111) . 
the observations cannot account for at least half of the predicted 
white dwarfs within the tidal radius, but observations and simu- 
lations agree for the number of young (<340 Myr) white dwarfs 
there. It is also worth to note that the simulations predict only a 
small number (five) of white dwarfs between 9 and 30 pc, which 
would be consistent with our findings of candidates in the 40 
pc range. On the other hand, the observations and simulations 
fully disagree for the older white dwarfs: the simulations sug- 
gest at least ten older (>340 Myr) white dwarfs within the tidal 
radius, but we do not observe any of them. It is very unlikely 
th at the ages of white dwarfs drawn fro m the cooling sequences 
of lWood 119951) : iFontaine et all d2001l) are significantly under- 
estimated given the consistency between white dwarf cooling 
ages and main sequence ages in the Hyade s and other clusters 
dDe Gennaro et all 12009; Jeff erv et afll2011|). We have verified 
that even with the calculations of Renedo et aD d2010), which 
are done by evolving white dwarf progenitors directly from the 
zero-age main sequence, the ages of the white dwarfs discussed 
in this work are only changed by ~10%. There are different other 



explanations for this discrepancy and we review some of them in 
turn. 

One possibility is that more white dwarfs than expected were 
kicked out of the cluster, and that most of them are now far 
from the center. T he kick velocities added to white dwarfs in the 
Ernst et all d201 ll) simulations can be seen as some kind of free 
parameter. However, since objects within the tidal radius ther- 
mally relax on short time scales, they can not be used to constrain 
this parameter. The observed spatial density and residual veloci- 
ties of white dwarfs outside of the tidal radius could in principle 
be used, but these numbers are poorly constrained by the cur- 
rent work. We suggest that new simulations could be run with a 
wide number of assumptions on the kick velocities, in order to 
predict ho w the number densit y and residual velocities would be 
modified. Fell hauer et al.l d2003l) suggest that older white dwarfs 
are preferentially lost due to their higher recoil speed, which is 
in turn due to the larger amount of mass loss of their high-mass 
progenitors, although it only remains an hypothesis. We note that 
similar simulations for the white dwarfs in other open clusters 
could be very helpful to constrain the distribution of kick veloc- 
ities. 

Another scenario is that fewer white dwarfs than predicted 
were ever formed in the Hyades, which would correspond to a 
change of the Initia l Mass Function (IMF) at its high-mass end. 
Dob bie et alj ((2006) estimate from their initial vs. final mass re- 
lation that the classical white dwarfs come from main-sequence 
progenitors in the range ~2.5-4 M Q . Therefore, to obtain more 
promising predictions from the simulations, one would have to 
modify the formation rate of white dwarfs for main-sequence 
stars with M > 4 M . Related to this, the fact that initial vs. 
final mass relati on for white dwarfs is not so well constrained 
for high masses (Dobbie et al. 2006, see their Fig. 3) might also 
have an effect on the number of white dwarfs vs. neutron stars 
and the age distribution of the remnants. Having a steeper initial 
mass function in this regime coul d have a strong effec t on the 
predicted number of white dwarfs ( B6 hm-Vitensdl2007l) . In this 
context we note that the present-day mass function of the Hyades 
derived in Paper I (Fig. 10) tolerates a steeper slope leading to 
less white dwarf progenitors. 

Thirdly, the topic of white dwarfs in binary sy s tems h as only 
been touched marginally in this paper. Holberg d2009h claims 
that 32 + 8 % of the local white dwarfs are in binary systems. 
Among the classical members, we find that three out of ten are 
binaries, a number similar to field objects. However, it would be 
dangerous to extrapolate that the fracti on of binaries is the same 
for field white dwarfs and the Hyades. Solli ma et alJ (|2010) esti- 
mate that the binary fraction of stars in open clus ters is between 
35 and 70%. Based on the Hi pparcos data from Perr vman et al] 
dl998l) . iBohm-Vitensel d2007l) found that the binary fraction in 
the Hyades increases from 30% for late-G to 88% for early- 
A stars. While we claimed previously that virtually all single 
white dwarfs have been discovered in the cluster center area, it 
is likely not the case for white dwarfs in unresolved binaries. 
For instance, V471 Tau and HD 27483 (Nos. 1 and 4) have 
not been discovered by systematic surveys, but from specific 
studies of Hyades main-s equence members dNelson & Youngl 
1970; Bo hm-Vi tense! 19931) . Since old and massive Hyades white 
dwarfs (e.g. with masses of ~0.8 M Q ) are completely hidden in 
the visible by their main-sequence companion in the majority 
of cases, they can only be detected by UV observations o f each 
main-sequence member. However, Bohm -Vitensel d 19951) found 
no white dwarf companion in a sample of IUE spectroscopic 
observations for 27 Hyades F stars. Furthermore. iPaulson et al] 
(2004) found no evidence of white dwarf companions in a radial 
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velocity survey of 94 Hyades stars. All in all, there is so far no 
observational evidence that the missing white dwarfs older than 
~340 Myr are hidden in binary systems. 

7. Summary 

The starting point of our work has been the Hyades white 
dwarf candidates kinema tically and photometrically selected by 
ISchilbach & Roseri d2012l) from the PPMXL survey. We evalu- 
ated the membership status of these candidates by using spec- 
troscopic and photometric observations, from which we deter- 
mined new distances of the white dwarfs from the Sun and from 
the cluster center as well as residual velocities with respect to 
the Hyades bulk motion. We find a very good agreement be- 
tween kinematic and spectroscopic distance predictions for the 
10 classical white dwarf members, confirming that the kinematic 
and photometric data sets, as well as model atmospheres that we 
rely on, are consistent with each others. 

We ruled out 6 of the Hyades candidates from Paper II as 
background stars, but also added 3 new objects by increasing the 
residual velocity upper limit for membership selection to 8 km 
s . There remain 14 candidates, five of which have cooling ages 
higher than the age of the cluster. We do not consider them as of 
Hyades origin. For another three, the atmospheric parameters are 
rather imprecise or could not be determined, leaving us with six 
spectroscopically confirmed hot and young white dwarfs consis- 
tent with a Hyades membership. In order to estimate the contam- 
ination fraction from field white dwarfs, we performed Monte 
Carlo simulations. The result is that ~2-3 hot and young white 
dwarfs may contaminate our sample yielding that the probability 
that all six are contaminants is low. A final confirmation of the 
candidates needs accurate radial velocity measurements properly 
corrected for gravitational redshift. 

The major discrepancy between observations of white 
dwarfs of Hyades origin and the predictions by simulations is the 
lack of about a dozen white dwarfs with cooling ages between 
340 and -600 Myr, which should be found within the present- 
day tidal radius (~9 pc) of the cluster. Possible explanations for 
this are: a) at their formation, within the first 300 Myr of the 
cluster's life, white dwarfs get a large (>5 km s~') kick velocity 
relative to the cluster motion which sweeps them farther out than 
40pc from the centre, b) the high-mass end of the IMF is steeper 
than the Salpeter slope of 2.3 assumed in the simulations, hence 
leading to a lesser number of massive stars, and c) there is still a 
number of so-far unrevealed white dwarfs hidden in binary sys- 
tems. There seems to be no way, at present, to decide which of 
these scen arios is the domi nant one in consistence with the con- 
clusion of lWilliamsl d2004l) . 
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